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Abstract. We consider the limit distributions of open quantum random walks on one-dimensional 
lattice space. We introduce a dual process to the original quantum walk process, which is quite 
similar to the relation of Schrodinger-Heisenberg representation in quantum mechanics. By this, 
we can compute the distribution of the open quantum random walks concretely for many examples 
and thereby we can also obtain the limit distributions of them. In particular, it is possible to get 
rid of the initial state when we consider the evolution of the walk, it appears only in the last step 
of the computation. 

III! 

1 Introduction 

Recently Attal et al. [21 S] introduced and investigated the open quantum random walk 
(OQRW) on graphs, which shows various different dynamical behaviors comparing to usual 
quantum walks and it includes the classical random walk as a special case. Then they 
considered the limit theorems for OQRW's: they have shown the central limit theorem [2]. 
The purpose of this paper is to further investigate the distribution and the limit theorem of 
OQRW's. 

A quantum analog of the classical random walk, called quantum walk (QW), has been 
intensively studied for the last decade (see [6l U\ IHl US]). The most remarkable difference 
between classical random walk and quantum walk appears in the central limit theorem. In 
the classical random walk, the limit distribution is Gaussian with scaling speed ^/n. On the 
other hand, in quantum walk the speed is linear in time, i.e., n, and moreover, the limit 
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distribution is far from Gaussian [H [5l El [9l [TO] . The main reason that makes quantum walk 
different from classical walk is the interference. 

The OQRW is different from the usual discrete-time QW. It was introduced in order to 
model the quantum efficiency in biological systems and quantum computing and it is based 
on the non-unitary dynamics induced by the local environments [31 H]. These random walks 
deal with density matrices instead of pure states. In [4j, Attal et al. developed the quantum 
trajectory approach for OQRW and by using this concept, they have shown the central limit 
theorem for OQRW's on the d- dimensional integer space Z*^ 

In this paper we focus on OQRWs on Z. We will introduce a concept of dual process to the 
original OQRW. By this we may think that the OQRW is a process where the environment 
evolves rather than the walker itself evolves as time goes on. This is particularly useful when 
we compute the distribution of the walk because the initial states are not relevant during the 
evolution. See section 2 for the details. Moreover, for many examples, we can compute the 
distribution of the walk very concretely and thereby we can also get the limit distributions of 
them. This paper is organized as follows. In section 2, we introduce the concept of OQRW 
following and new concept of dual process. Next we give our main result. Section 3 is 
devoted to the proof of our main result and some preparation that is useful for examples. In 
section 4, we consider several concrete examples. 



2 Open Quantum Random Walks 

In this section, we give a brief definition of OQRW and define a dual process of it. Then we 
state the main result. In order to compare, we first shortly review the usual quantum walks, 
so called unitary quantum walks. 



2.1 Unitary Quantum Walks 

The discrete-time QW is a quantum version of the classical random walk with an additional 
degree of freedom called chirality. The chirality takes values left and right, and it determines 
the direction of the motion of the walker. At each time step, if the walker has the left 
chirality, it moves one step to the left, and if it has the right chirality, it moves one step to 
the right. In this paper, we put 



1^) 



1 




\R) 




1 



where L and R refer to the left and right chirality state, respectively. 

For the general setting, the time evolution of the walk is determined by a 2 x 2 unitary 
matrix, f/, where 



U 



a b 
c d 



with a, 6, c, (i G C and C is the set of complex numbers. The matrix U rotates the chirality 
before the displacement, which defines the dynamics of the walk. To describe the evolution 
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of our model, we divide U into two matrices: 



a b 




Q 




c d 



with U — P + Q. The important point is that P (resp. Q) makes the walker moves to 
the left (rcsp. right) at position x at each time step. For example, the Hadamard walk is 
determined by the Hadamard gate U = H: 



H 



1 

7^ 



1 1 
1 -1 



The walk is intensively investigated in the study of the QW. 

Let S„(Z,m) denote the sum of all paths starting from the origin in the trajectory con- 
sisting of I steps left and m steps right at time n with n — I + m. For example, 

S2(i, 1) = QP + pg, 

S4(2, 2) = Q2p2 + P2Q2 ^ QpQp ^ pgpg ^ pg2p ^ gp2g 

Let Z+ = {0, 1, 2, . . .}. The probability that our quantum walker is in position a; G Z at time 
n e Z+ starting from the origin with ip = *[q;, (3] with a, /3 e C and \a\'^ + = 1 is defined 

by 

P{Xn = x) = ||E;„(i,m) (fW^, 



(1) 



with n — I -\-m and x = — Z + m where /3] means the transpose of [a, p]. 

The reason that we call this quantum walk the unitary quantum walk is that the evolution 

of the walk is the unitary transform of a state in a Hilbert space. To say more concretely, 
let us denote by "He := the space of intrinsic structure, namely the chirality and let 
ICs '■= ^^(Z) the space of positions. The Hilbert space on which the quantum walk evolves 
is given by 

and any state, i.e., a unit vector of T-L is given by 

Let us denote by T the left translation in Z^(Z): 

{Ta){x) = a{x + 1), for a = {a{x))xez- 
T is a unitary map whose adjoint is the right translation: 

{T*a){x) = a{x - 1), for a = {a{x))x&z- 

Let ^0 £ ^ be any initial state. Then the dynamics of 1-dimcnsional quantum walk driven 
by the unitary operator U — P + Q in the above is represented as 
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where ipn = ('?/'n(a;))xez G "H is the state at time n. It is easy to see that the operator 
PT + QT* is a unitary operator (here tacitly we understand the operators P and Q are 
the natural extensions of the original ones to ®x&^iHc) and hence the quantum walk is a 
unitary evolution on the Hilbert space 1-L. The probability density of (|5]) is then nothing but 

P{Xn = x) = Un{x)\\\ 

with initial state V'o = ® |0). Here {|a;) : x G Z} denotes the canonical orthonormal basis 
of /2(Z). 

The limit distribution of quantum walk is very different from that of classical random 
walk. It is ballistic instead of diffusive in the sense that Xn/n converges weakly to a limit 
whose density function, rigorously shown by Konno [Sl[in], is given by (in the case ahcd ^ 0) 



7r(l — x^jy |ap — 

where /3 is a constant depending on the initial state if and the unitary matrix U . Here 
Ia{x) = 1 (x e A), = (x ^ A). 



2.2 Open Quantum Random Walks 

In this subsection we introduce OQRW's following [3]. The OQRW's can be defined on 
any dimensional integer spaces as well as on any graphs, but here we confine ourselves to 
1-dimensional space Z for simplicity. 

As before He = is the Hilbert space for the intrinsic structure and ^5 = /^(Z) for 
positions. Let B and C be two linear operators on Tic, i.e., 2x2 matrices, such that 

B*B + C*C = I. 

Define a completely positive map on the density matrices of Tic by 

C{p) := BpB* + CpC*. 

The OQRW lifts this map to "H = Tic ® Tis in the following way. We consider density 
matrices on "H of the type: 

P = X^Px- ® |a;)(x|, 
where each is a positive matrix and satisfy 

J2T^r{px) = ^- 

Let C be an operator that maps on such density matrices as follows: 

C{p) = J2iBp,+iB* + ® |x)(x|. (2) 
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The OQRW is an evolution obtained by iteration of the map £. Let p(o) :=po® |0)(0| be 
the initial state, i.e., a density matrix on %. Then, the evolution of OQRW on Z generated 
by B and C is given by 

X 

pi"+^) = Bpi^lB* + Cpi%C\ xeZ, n = 0,1,2,- ■■. (4) 

The probability distribution to find out the walker at site x at time n is given by 

pi")=Tr(p^")), xGZ, n>0. (5) 

In Attal et al. introduced the concept of quantum trajectory to OQRW's, from 
which they could show the central limit theorem for OQRW's f2]. Applied to OQRW's on 
1-dimensional space Z, it says that one can introduce a Markov process on 8{'Hc) ^ ^5 
where SiTic) is the space of density matrices on Tic such that the distribution of the 
space component of the process coincides with that of OQRW. To say more in detail, it is 
a Markov chain [pn, Xn)n&N with values in Eijic) x ^ with the following transition rule. 
From any position (p, X) it jumps to one of two states: to {—BpB*,X — 1) with probability 
Pb '■= Tr^BpB*) or to {^CpC*,X + 1) with probability pc '■= Tr(CpC*). Then the central 
limit theorem shown in [2] reads as follows (stated in 1-dimensional case). 

Theorem 2.1 ([2], Theorem 5.2) Consider the stationary open quantum random walk 
on Z associated to the operators {B, C}. We assume that the completely positive map 

C{p) = BpB* + CpC* (6) 

admits a unique invariant state poo- Let {pn, Xn)n>o be the quantum trajectory process to 
this open quantum random walk, then 

Xn — nm 

converges in law to the Gaussian distribution A^(0,(T^) in R, with mean 

m = TiiCp^C*) - TiiBp^B*) 

and variance 

= Tr{Bp^B* + Cp^C*) - + 2Tr[{Cp^C* - Bp^B*)L] - 2mTr{p^L), 
where L is the solution of the equation 

L - £*{L) = C*C - B*B - ml. 
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2.3 Dual Processes and Main Result 



In this subsection we introduce the concept of dual process to the OQRW and state our main 
result. Recall that 1-ic = and so B{l-Lc) = M.2, the algebra of all 2 x 2 matrices. From 
now on we regard A^2 as a Hilbert space equipped with the Hilbert-Schmidt inner product: 

(v4,5) := Tr(A*5), A.BeMi- (7) 

Let Ai := ®x&M.2 be the direct sum Hilbert space. Recall the left and right translation 
operators T and T* on The operators T and T* naturally extend to Ai. We will 

use the same symbols whenever there is no danger of confusion. Let K := (— vr, it] and we 
understand it as a unit circle on the plane. The Fourier and inverse Fourier transforms 
between /^(Z) and L'^{K, ^dx) are defined as usually: 

a{k) := ^e-^'='^a(x), a = (a(x))^6z G /^(Z), 

/» ■■= ^ f e''''f{k)dk, feL\K,^dk). 
Ztx Jj^ Ztt 

For each k E K, let ^Ak be the copy of the Hilbert space A^2 and let 

/•® 1 
Jk 27r 

be the direct integral Hilbert space. The Fourier transform also naturally extends to the 
transform between A4 and Ml: for A = {A{x))x<^z ^ -M, 

A ■= {A{k))keKeM, 
A{k) := ^e-^'^^A(x), 

and similarly for the inverse transform. 

Now let us introduce the left and right multiplication operators on the Hilbert space A^2- 
For any i? G A^2, the left multiplication Lb and the right multiplication Rb are defined on 
M2 by 

Lb{A) := BA, Rb{A) := AB, A e M2. (8) 

Notice that (Lb)* = Lb* and (Rb)* = Rb*, and for any B and C in A^2, Lb and Rc 
commute: LbRc = Rc'Lb- However, Lb and Lc do not commute in general. The operators 
Lb and Rb are positive definite if 5 > 0. Without mentioning further, we will use the same 
symbols Lb and Rb for the extensions to Ai. Thus, as an example, for A = {A{x))xez ^ -M., 

LbRc{A) = {BA{x)C)x^z G M. 

Let us now come back to the OQRW's on Z. The state at time n can be understood 
as an element of Ai. Then the dynamics of OQRW's becomes an evolution on Ai given 
by 

p("+i) = {LbRb'T + LcRc'T*)p'^^\ (9) 



6 



Therefore the solution to (Q becomes simply 

p(") = {LbRb'T + LcRc^TTp^'^- (10) 
If we look at the evolution in the Fourier transform space, then it becomes 

^){k) = [e'^^LBRB'+e-'^LcRc^T^^k). (11) 

Notice that p^^\k) is the constant (operator valued) function po because p^^^ = Po ® |0)(0| 
(of course we can take quite general initial state p^^^ not localized at the origin). 
In order to get the probability distribution, let us define a "dual process": 

Definition 2.2 The dual process to the OQRW generated by B and C is the process 
Yn := {Ynik))keK e M defined by 

Yn{k) := {c'^Lb^Rb + e-'''Lc*RcT (^)- (12) 
The reasoning for the nomenclature becomes clear from the following relation, which is our 
main result. 

Theorem 2.3 The probability distribution of the OQRW is given by 

Pi"^ = ^ X ^''''Tr (poYM) dk. (13) 

Remark 2.4 Notice that in the formula for the distribution of the walker in the above the- 
orem, the initial state po does not change at all as time goes on. Instead, the environment 
denoted by B and C evolve. 

The easy proof of Theorem 12.31 will be given in the next section. 



3 Proof and Some Analytic Preparation 

In this section we provide with the proof of Theorem 12 . 3 1 and we give some analytic prepara- 
tion which will be useful when we consider asymptotic behavior of some functions. We start 
with 

Proof :[o{ Theorem 12.3] Recall the inner product in ([7]). By using Fourier transform and the 
formula (fTTj) . we see that 

= Tr(p(")) 

= (/2,pi")> 

- ^Xe-(/„p^(.))ci. 
1 

2^ 



]- [ e'''- {{c-'^Lb^Rb + e'^Lc^RcT ih), Po) dk 
Jk 

]_ I e'k.rj.^{po{e'''LB*RB + e-''Lc^RcT{l2))dk 
Jk 



27r 

1 ' Akxr 



— / e''^Ti{p^Yn{k))dk. 



K 



□ 



In the next section we will consider several examples and compute the distribution concretely 
by using Theorem 12.31 

Now we introduce some analytic result which is not only interesting in itself but also 
useful for studying an asymptotic behavior of functions. But, it may be well known in 
analysis. 

Proposition 3.1 Let [a,b] C M &e a finite interval and let / : [a, 6] — ?■ M 6e a continuous 
function such that \ f\ has a unique maximum at a point c G [a, b]. Then for any continuous 
9 ■■ [a, b] -> R, 

1 '■^ 



lim — / f{xYg{x)dx = g{c), 

ra-S>oo Ctr 



where an = f{x)'^dx. 



Proof: By multiplying a constant, if necessary, we may assume maxx£[a,b] \f{x)\ = /(c) = 1- 
Also, we may assume g{x) > 0. We will assume c G (a, 6), because the proof for the case 
c = a or 6 is similar. Given a sufficiently small e > 0, we can find a 5 > such that if 
|x — c| < 6, then f{x) > and \g{x) — g{c)\ < e. Moreover, if — c| > 6, then < 1 — e. 

Also, we can find a < 6i < 6 such that if |x — c| < 6i, then f{x) > 1 — e/2. First we show 
that there is an no G N such that 

a„ > if n > uq. 

In fact, we can divide the integral 

an= [ f{xYdx= I f{xYdx+ I f{xYdx. (14) 



Notice that 



/ 



f{xfdx > [ f{xYdx > (1 - e/2)"25i 

5} J{\x-c\<5^} 

|/(x)|"cix < (l-e)"(6-a). 



'{\x-c\>5} 

Thus we see that 



By and ([U]) we have 



^.^ J{\x-c\>s} \J y n ^ Q ^^g^^ 



lim—/ /(x)"cia; = 1 and lim — / \f{x)['dx = Q, (16) 

an J{\x-c\<5} "n J {\x-c\>5} 

and the claim follows. Now notice that 

— / f{xYg{x)dx = — I f{xYg{x)dx + — f f{x)'''g{x)dx. 

Ja Oin J {\x-c\<5} J {\x-c\>5} 
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Since \g{x)\ < M for some M > 0, for n > riQ 
^ ' f{xYg{x)dx 



J{\x~c\>6} 

Hence by f lTB|l we see that 



< M— [ \fix)\''dx. 

O^n J{\x-c\>5} 



lim — / f{xYg{x)dx = 0. 

n^oo an J{\x-c\>5} 

On the other hand, for n > no we have 

^^^^ / fixTdx < — [ f{xrg{x)dx < M±I 

O^n J{\x~c\<S} Oin J {\x-c\<5} '^n J {\x-c\<5} 



< 2e. 



(17) 



f{xydx. 



Therefore, by f[T6l) we see that for sufficiently large ra's 



1 



'^n J {\x-c\<S} 

By (|T7|) and f|T8|) . the proof is completed. □ 



f{xfg{x)dx - g{c) 



4 Examples 

In this section, in order to see the usefulness of our theorem we consider several examples. 
We also compare with the result of [2]. 



4.1 Example 1 

For p and q such that p + g = 1, p, g G [0, 1], let 



' 1 " 




' " 


^_ 


c = 






B 



Since S and C are all diagonal matrices. Lb and Lc, -R_b and i^c commute. So, it is very 
easy to compute 



Yn{k) = {d^LB^RB + e'^'^Lc^RcYih) 

n 



1=0 
Akn 



1 




ik(n-2l) 


'1 


0" 




'0 


0" 







p. 











1 
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Therefore, since po 



a 
b 



, by Theorem 12.31 we have 

^TT Jk 

Thus a standard argument imphes 

Proposition 4.1 Consider OQRW with initial state po- //p ^ (0, 1), then as n —> oo, 



n 



a + 6 (5, 



9-P' 



(U) 



Xn-{q- p)n 



n 



bN{0,4pq) 



We notice that in this example the equation (|6]) has infinitely many invariant states. 

4.2 Example 2 

We take 



B 



' bn 


" 


, c = 


' 


Cl2 


_ b2i 








C22 _ 



(19) 



such that U = B + C G U{2), where U{2) is the set of 2 x 2 unitary matrices. We suppose 
that 



Let Pi : = 



1 


easily shown that 



P = \bu 
and P2 : 



P = l-|&2lP 



I |2 _ 1 I |2 
C22 — i — C12 . 




1 



be the projections. From the definition of B and C it is 



and 



B*B = Pi, B*PiB = pPi, B*P2B = (1 -p)Pi 

C*C = P2, C*PiC = (1 - p)P2, C*P2C = pP2. 



(20) 
(21) 



Therefore, Yn{k) = {e^^LB*RB + e ^''Lc*Rc)'^ (h) becomes a linear combination of Pi and 
P2. 
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Lemma 4.2 Let 

Then the coefficients are determined by 



'ik I 



l-p) 



l-p) 



Proof: We have Yn+i{k) = [c'-^Lb'-Rb + e ^^Lc*Rc) yn{k). Inserting 
relations (120|) and (l2T|l we get the recurrence relation 



(n+l) 

n 

(n+2) 



(k) 



l-p) 



e ^^p 



p) 



'af\k) 
a^^\k) 



Since 



'af\k) 
af\k) 



the result follows. □ 



(22) 



and by using the 
(23) 



Proposition 4.3 The OQRW generated by B and C in fjlQl) is a correlated random walk 
and as n ^ oo, 







p 



' 1 — p/ 



Proof: By Theorem 12.31 and Lemma 14.21 we have 

p^^^ = ap^^\x) + hpl^\x), 

where 



(24) 



pf\x):=^ / e^^^af\k)dk, j = l,2. 



Since 



-ik I 



e p e (1 — p) 
l-p) 



1 

^ik Q 

e 



-ik 



-^-ik Q - 




aS"+'Hfc)" 




e'^ 









p 1 — p 

1 — p p 

p l-p 

1 — p p 



from (123|) we get 

aS")(fc)' 
af\k) 



We multiply e^^^ to the both sides of the above equation and integrate w.r.t. k. Then we 
obtain 



p'f''^^\x — 1) 
pt^'\x + l) 



p l-p 
l-p p 



Pi [x) 

p\ '{x) 



or rewriting it we get 

p("+i)(a;) = ppf^{x^l)^{l-p)p^^^{x^l) 

P2"'^^'*(x) = (1 —p)p^^^{x — 1) ^pp^2^{x — 1). 

From the above relations and fl2^ . we get the result (see [12], for example). □ 

We see that this example falls into the class that the result of [2] may be used. We 



compute that (E]) has a unique invariant state pc 



1 (&11&21 + C11C22) 

,(&11&21 + C22C12) 1 
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4.3 Example 3 

Let us define 

B = 



1 



\ P-\ 



c 



7 



xh-\ 



with p + q = 1, p,q E [0,1] and < 7 < mm{^/2p, y/2q}. It is promptly shown that 

B*B = Pi+pP2, B*PiB = Pi, B*P2B = pP2 (25) 

and 

C*C = qP2, C*PiC = 7'P2, C*P2C = qP2, (26) 

where p = p—\ and q = q — \- Therefore, F„(A;) = (^e^^LB*RB + e~^''Lc*Rc)^ (h) is again 
a linear combination of Pi and P2. By denoting F„(A;) = af'\k)Pi + a^^\k)P2, we repeat 
the method done in Lemma [4.21 using the relations fl2^ and fl2B]) . Then we easily get 







e *^7^ e*^p + e 



_N I 



(27) 



Lemma 4.4 VKe have 



P: 



n—1 j 



j=0 1=0 



n 



j=0 



Proof: Recall that Yn{k) = a^^'\k)Pi + a^^'>{k)P2 with ^(A;), j = 1, 2, being defined by 
By Theorem 12.31 we see that 



An) 



(«)/ 



1 



p-'=a— / e''="a!"^(A;)rfA; + 6-^ / e'^'=a^2\k)dk. 



1 



(28) 



Inserting the formulas for a^"'\k), j = 1, 2, in (^7^ into f l2Sl) . we easily get the result. □ 



Proposition 4.5 Consider OQRW with initial state po 



a 
6 



® |0) (0|. r/ien 



n 



5-1 



as n ^ 00. 
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Proof: Let us compute the characteristic function 0x„/ra(^) '■= E (^e**"^"/"). By using Lemma 
lOl we get 

as n — )• oo. In the last hne we have used the fact that (p + ge^**/")" — )■ because 
|p + ge^**/"| < p + g < 1. Also similarly (pe"^**/" + g) " -> 0. Notice that e~** is the 
characteristic function for the distribution □ 



Remark 4.6 By the same method as above we can in general show that for all a > as 
n — !■ oo 

If we rely on Theorem \2.1\ we can show that m = —1 and = and the result says that as 
n ^ oo 



with cr^ = 0. 



4.4 Example 4 



Here we will consider an example of OQRW whose distribution is a mixture of normal 
distributions. Let < e be a small number such that 2ea(e) < 1/2, where a(e) := a/1/2 — e^, 
and let 6* G M. Define 



B 



' a(e) ee^^ " 


c = 


ee^^ a(e) 





a\t) —ee 
-te'^ aie) 



(29) 



It is straightforward to see that all the matrices 5, 5*, C, and C* commute with each other 
and 



B*B = 
and thus 



1/2 2ea(e)cos^ 
2ea(e)cos0 1/2 



C*C 



1/2 -2ea(e)cos^ 
-2ea(e)cos^ 1/2 



(30) 



B*B + C*C = h. 
The two matrices B*B and C*C are spontaneously diagonalized as 



B*B = U 
where the eigenvalues are 



A+(e,^) 
A_(e,i 



U and C*C = U* 



A_(e, 






A+(e,i 



X±{e,9) = l/2±2ea(e) cos^ 



(31) 
(32) 
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and U is the unitary matrix given by 



V2 



1 1 
1 -1 



(33) 



Because of the commuting properties of the matrices, we can easily compute that 

Y^k) = {e'^LB*RB + e-'''Lc*RcT{h) 
= {e'''B*B + e-'^C*CY 



E 

1=0 
n 

E 

1=0 



n 



n 



e-'^^''-'^^\B*B)\C*C) 



-ik{n—2l) JJ-* 



n—l 



u. 



(34) 



Now by using Theorem 12.3^ we get the following result. 



Proposition 4.7 Let B and C be the matrices in fl5CT]) . X± = X±{e,6) the eigenvalues of 
B*B and C*C in fl52]) . and let U be the unitary matrix in fl5^ . Let po be the initial state 
on M.2- Then the probability distribution of the OQRW defined by B and C is given by 

in) M I (n) 

Px ' = +«2P2,x, 

where ai = {UpoU*)u and a2 = {UpoU*)22, and p^^l, j = 1,2, are the distributions of 
random variables Xj^n, j = 1, 2, respectively, whose asymptotic behavior are as follows: as 
n — )■ oo. 



Xi^n - (A- - \+)n 



n 



X2,n - (A+ - A_)n 



n 



iV(0,4A+A.), 
A^(0,4A+A_). 



Proof: From we see that 

TiipoY^ik)) = J2 ('^)e-'=("-2')(aiA^Ar' + a^ALA^O, 
1=0 ^ ^ 

where Oi = {UpoU*)ii and a2 = {U PqU*)22- Therefore by Theorem 12. 3[ we have 

i—n V / 



Now the result follows from the standard arguments. □ 
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4.5 Example 5 

Define 



B 



1 

7! 



1 1 

1 



c 



1 

7! 



1 

-1 1 



We let 



Yn{k) 



(e^'^LB^RB + e-^'Lc^RcTih) 



By directly computing we get the recursion relation: 



a 



"11 

(n+l) 

at'\k) 



(k) 
(k) 



\k) 



Y{k) 



4? (A:) 



where 



2 cos A; 



— e 

2 COS A; 




— e 


2 cos k 



-ik 
,—ik 



-ik 



—e 
2 COS k 



Thus we have the solution 







"1" 


atl\k) 


= r(A;)" 










1 



By Theorem 12.31 we have 



^ / e*'=^Tr(poV;(A;))rfA; 



^ J e'^"" (aa^^^k) + 64?(A;)) dk. 



27r 

We can introduce a combinatoric way to compute the distribution. Notice that 



1 n 
n + 1 



and C*"C^ 



+ 1 
— n 



-n 



Thus 



+ 2 



/9. 



We notice also that 



TrC*"C" 
15 



Since Yn{k) = {e^^LB*RB + e ^^Lc*Rc)'^ (h), by expanding the power and using Theorem 
12.31 we see that pi""* is the sum of all the contributions from the terms of the type 

{Lc.RcY'iLB^RBr'iLc^RcY'-'iLB^RBY'-' ■ ■ ■ {Lc*Rc)''{LB*RBT{h) 

where Yl]=ii.h ~ ''^j) ~ ^ '"i — — 1' i = 2, ■ ■ ■ , s, Z-,- > 1, j = 1, ■ ■ ■ , s — 1, 

and Is > 0. We would like to expand (j^. By (jH]), the term B*''^B'''^ in the middle is 
equal to -^r-^'^ '~ C*'^^C'^'^. Inserting this into fj42|) . we get a sum of two sequences, whose 
middle terms are C*''^C^^ multiphed by a factor and (7*('i+''i)(7('i+''i) multiphed by a 
factor —1, respectively. We continue this process successively. For it, it is very convenient 
to understand (1421) as a random walk path (assume all r^'s and Z^'s are greater than for 
simplicity): the walker goes upward ri units, then goes li units downward. Then it goes r2 
units upward and I2 units downward, and so on. So, the path is a continuously connected 
lines consisting of 2s segments (which have different lengths of ri, /i, etc.). We will further 
simplify the notation by denoting it just as a sequence (Z^, r^, ■ ■ ■ , Zi, ri). We will make short 
the sequence step by step by applying the process of "cutting" or "unfolding". For example, 
at first step, if we make cutting we will get the sequence (Z^, r^, ■ ■ ■ , Zi) with a weight 
Instead, if we make unfolding at the first step, we get a new sequence {ls,rs, ■ ■ ■ ,h + ri) 
with weight —1. Any operation shortens the sequence by length 1. We continue this process 
until we get a length one sequence, or for the random walk path, until it remains a single 
segment of length Z^ + Z', say. The resultant matrix is nothing but C**''"*"' (jh+v g^j^^^ need 
to compute the trace Tr(poC*'"+''C'^+''), which is simply (a ((Z^ + I'Y + 1) + 6). We 

summarize this process as a theorem. Below Yl^h--- h r-i ••• r-s i^^ans the sum over sequences 
such that Yli]=i{h ~ ^j) ~ ^ > 0, Tj > 1, j = 2, ■ ■ ■ , s, Ij > 1, j = 1, ■ ■ ■ , s — 1, and 

Is > 0. CU{ls,rs, ■ ■ ■ ,h,ri) means the set of all sequences of shortening process of cutting 
and unfolding upto a single term and for tt G CU{ls, ^s, ■ ■ ■ , h, ri), Z(7r) is the length of the 
remaining single segment for the process vr, and cj(vr) is the product of the weights of vr 

Theorem 4.8 The probability distribution for the OQRW determined by B and C in (l35|) 
is given as follows: 

'i,--- ,ls,ri,--- ,rs TT&CU{ls,rs,--- ,li,ri) 

where C = Cifls^QandC = Bifls = Q. 

As an example, let us compute the distribution of X4 for the case a = 6 = 1/2. Since 
Ti{B*^B^) = Tr(C*'C') for all Z > 0, it is easy to see that the distribution under the 
assumption is symmetric. So, we only need to compute P(X4 = 4) and Pi^X^ = 2). The 
random walk path leading to X4 = 4 is a single segment consisting of 4 upward units. Or, 
in the symbol of finite sequence, it is just (ri) = (4). Thus, we get 

P{X, = 4) = Tt{p,B*'B') = = 1. 
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For X4 = 2, we have 4-paths: (- + ++), (H h+), (+ H h), (+ + H — ), or in symbols 

h ri r2 h ri r2 h ri h 

of sequences (1,3), (1,1, 2), (2,1,1), and (3,1), respectively. For each symbol we apply 
cutting-unfolding process. 

(1, 3) : ^^TripoC*C) - Tr(poC*^C^) = ^ 
(1,1,2): ^!±^(1,1)-(1,3) 

= (^^Tr(po5*5) - Tr(pofi*'fi')) - (^^-^MpoB* B) - T,{p,B*'B'] 

1 

~ 54 

(2,1,1) : (2,1) -(2, 2) 

= (Tr(po5*'fi') - Tr(po5*'fi')) - {^^^Tiip.B*' B') - Tr(po5*^i?^)) = ^ 

(3, 1) : TiipoB^'B') - Tr(po5*'i?^) = 

54 

Thus summing all the contributions we get P{Xi = 2) = |. Using the symmetry, we see 
that p4, the distribution of X4, is equal to 

1 2^ 3^ 2^ 1 

^4 = 9^-4 + g<^-2 + g^o + + -54. 

From now on we discuss the asymptotic behavior of the distribution. Recall the matrix 
Ynik) in (15^ and its representation in (155]) . The eigen-equation of Y{k) is 

2cosfc\ /,o ,,n 4cos^/i;+l 2 cos A;(4 cos^ A; + 5) \ 
A - — j (a^ - 2 00s *A^ + 4^ ^ j = 0. 

SO the eigenvalues of Y{k) are 

2 cos , 2 cos k 1 f ^ 1\ , ^ . 

A, = ^+i|(-l+,V3)^ + i±i^i. A3 = A.. (44) 



where 



Let 



^ = ^(A;) = (^2cosA; + V4cos2/c + ly^^ (45) 



— Ai — Aq, A2 — X2 — Aq, A3 — A3 — Aq. 
Then a direct computation gives 

Y{k) = S-^ diag[Ao,Ai,A2,A3] S, 
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where S := with 



R 



1 
6 



2 - ee'^'Ai 

6 1-9AI 

-6 1 - 9Al 

-2 + 6e-*^Ai 



2 - 6e^*^A2 
1-9^2 



1 - 9Al 
-2 + 6e-^^A2 -2 + Ge-^'^Ag 



2 - Ge^'^As 
1-9^2 
1-9AI 



and diag[ao, ai, . . . , a„] denotes the diagonal matrix whose {i, i)-component is Oj. Here 

R-' = 

W31 W32 W33 W34 

_ W41 W42 W43 m;44 

Recall the density formula in fl39|l . We need to compute a^Jj'{k), k = 1,2, which we can 
obtain from the diagonalization of Y{k). We note that deti? = 4a/3(4 cos^ k + 1) sinfc/9. 
After a little computation we have 






1/2 


-1/2 





W21 


W22 


W23 


W24 




W32 


W33 


W34 


W41 


W42 


W43 


W44 



1 

2^ 
1 

2^ 



Akx 



Akx 



aa^;i\k) + ba^^{k)jdk 
-2ismk Aj (wj+i^ia + Wj+i^J)) A" j dk. 



(46) 



In order to get an information of the asymptotic behavior or as n — )■ 00, we need to 
investigate the eigenvalues more carefully, and then we will rely on Proposition 13.11 For that 
purpose we will rewrite the eigenvalues. Recall the eigenvalues in (H3|) and (jB]), and the 
function ^(fc) in (H5|) . Since cos A; appears in the eigenvalues, we let m := cos fc. Then u varies 
in the interval [—1, 1] and we have 



Further, we define 



^ = ^{u) = (2n + V4u2 + 1)^/1 
s = s{u) := C,{u) r^^, — 1 < M < 1. 



It is not hard to show that 



-2 + v^)i/3 



V5 



(2 + V5y 



/3 



2 

75 + 1 



Moreover, on the interval [—1, 1], the function s{u) is increasing with 

s(-l) = -1 and s(l) = 1. 

We can also check that 



Au. 
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Therefore, the eigenvalues can be rewritten as 

Ao 
Ai 



6 



A2 

A3 

Since — 1 < s < 1, we see that 





a;. 



|Ao|<2/3, |A2|<v^- 



(47) 



Only the eigenvalue Ai moves fully on the interval [—1, 1]: Ai(m = —1) = —1, Ai(m = 1) = 1. 
Regarding K := (— vr, tt] as a unit circle in the plane, as usual, it is not hard to see that the 
eigenvalue Ai = Ai(A;) is anti-symmetric in the sense that Xi{k + tt) = — Ai(/c). Moreover, it 
behaves very much similar to the cosine function. In particular, Ai(A;) > on [— 7r/2,7r/2] 
and it is negative on K \ [— 7r/2, 7r/2]. Let us define a scaling constant an by 



Ctr. 



■k/2 



tt/2 



\i{kYdk. 



(48) 



Lemma 4.9 For j = 0,1,2,3, let gj{k) he continuous functions onK. Then 

lim — / gAk)\AkYdk = Q, 
ttn Jk "Zlf 

and 



lim / g^{k)\^{kf^dk = g^{Q)+g^{Ti), 

n^oo a2n Jk 

lim^— / g,{k)\,{kf^-^dk = g,{Q)-g,{'K). 



Proof: First, since Xi{k) is continuous and Ai(0) = 1, as in the proof of Proposition 13. H it 
is very easy to see that for any a/2/3 < g < 1, 



lim — = 0. 



n-!>oo ar. 



(49) 



On the other hand, by (j4 



K 



Y^g,{k)X,{krdk 



0((2/3)"/2). 
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From this and by ( H9|) . the first assertion follows. For the second assertion, we divide the 
integral: 

— ! gi{k)\i{kYdk = — I gi{k)\i{kYdk + — I gi{k)\i{kYdk. 

O^n Jk Oin J[-n/2,7T/2] Oin J K\[—k /2,-n /2] 

Noticing the anti-symmetry of Ai(fc), i.e., \i{k + tt) = — Ai(A;), the result follows from 
Proposition 13.11 □ 

Looking at the formula fl46l) . by Lemma [4.91 we see that asymptotically the term containing 
A" dominates. By direct computation we have 

^21 = f ^e-^' - (^2 + ^3) + 3e-*%A3 



deti? V3 
I2-A3) fl^ 
detR VS' 

Now we can get the proper asymptotics for the density 

Theorem 4.10 As n ^ 00, the asymptotic behavior of p'^^ is as follows. 

, Paf"'' I 1/^5 if X is even, 

lim = < 

n^oo I 0, if X is odd. 

p^x^~^^ )0, if X is even, 

lim = < 

n^oo a2n-i if X is odd. 

Proof: The proof will follow from Lemma 14.91 with 

gi{k) = C{k){aBr{k) + hB2{k)). 



where 



C{k) 



-9ie*'=^(A2 - A3) 
47rv/3(4 cos^ k + 1) 



Biik) = Ai Qe-*'^ - (A2 + A3) + Se-^MsAs 
B2ik) = Qe^^ - (A2 + A3) + 3e^%A3 

We need to know the values gi{0) and (71 (vr). Let us define a symbol rj by 

V{0):=1, r]{n):=-l. 

By directly computing, we get 

3 3 

C(0) = — and Cfvr) = — cosvrx. 
An An 
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Also it is easy to check that when restricted to {0,7r}, Ai = ^r] and the quantities in the 
brackets (■ ■ ■ ) of -Bi and B2 are equal to 2r]. Thus we see that when restricted to {0, vr} 

B^ = B2 = 2/3. 

Combining these we use Lemma [4.91 to get the result. □ 

Concerning the central limit theorem of this example we have the following result. 

Theorem 4.11 For the example of this subsection, we have as n 00 




N{0,8/9). 



Proof: Let us consider the characteristic function 

Here, pi""* is given in ( H^ . but since the eigenvalue Ai dominates we have 

/7r 
-TT 

where 
with 



By putting y = xj ^fn and taking a change of variable m = y^k we have 

Notice that the function g{k) is bounded and continuous and Ai(A;) behaves very much similar 
to cosine function on the interval [— 7r,7r]. Now we expand the interval to [—^Jn^T^ ^Jn^^\ 
and take a power n to Ai. As n grows, the function Ai(m/A/n)", when integrated with a 
multiplication by a mild function gimj ^fn\ will pick up the values of gimj ^fn) at m = 
and m = ±y/nTT with dominating factors of itself. Notice that 5^(0) = 5'(±vr) = We first 
consider the behavior at m = 0. For that we notice that Xi{k) has a Taylor expansion at 
/c = as 

x,ik) = i-h' + ^^k' + o{k'). 

Therefore, 
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Thus, the contribution to E[e**"^"/^] is evaluated as 

1 1 



oo 

4^2 



1 



oo 

1 ,,2 



-)■ ^ / e**^e "^'^"dy as n -> cx), (50) 
v27ro"^ J-oo 

where = 8/9. Next we consider the effect coming from the factor \i{k)^ aX k = ±7r. For 
this, it is convenient to shift the integration interval as 

Jo 

Now Ai(7r) = —1 and by a similar argument as above the contribution to E[e**"^"/^] is 

1 1 r°° 

= (-1)" V e^'y^—e'^'^y-y/^e-Tey' 
^-^ y/n 2iT 2 



Now as n — )■ oo, by an argument of Riemann-Lebesgue Lemma, the last term converges to 
0. Combining this with (150|) . we conclude that in the example 4, 

^" ^ Arm, 8/9). 
'n 



□ 

This example was also dealt with in [2] too. There they computed the invariant state of 
([6]) obtaining poo = They also computed the mean m = and variance cr^ = |, the same 
result as we obtained here. 
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